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Abstract. 

A systematic analysis of resonant x-ray Bragg diffraction data for UPda , with signal 
enhancement at the U M/y edge, including possible structural phase-transitions leads 
to a new determination of the space groups of the material in the phases between 
To = 7.8 K and T+i = 6.9 K as P222i and between T-i = 6.7 K and T2 = 4.4 K as 
space group P2i. In addition, the quadrupolar order-parameters, (Qab), inferred from 
diffraction data in the phase between T^i and T2, are {Qxz) and (Qyz) at the (103) 
Bragg reflection, and (Qxy) at the (104) reflection. 
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1. Introduction 

Phase transitions in materials are often driven by a cooperative action among degrees 
of freedom of the electrons. The paramagnetic-magnetic transition, for example, may 
arise from interactions between magnetic dipole moments. The potential importance of 
multipolar, e.g., quadrupolar, interactions in f-electrons system has long been recognized 
[H m [3]. Experimental techniques utilizing resonant x-ray Bragg diffraction (RXS), 
which have evolved in the past decade, are now the preferred techniques for exposing the 
interplay of charge, spin and orbital electron degrees of freedom in complex materials. In 
short, RXS reveals multipoles with a directness of purpose not available with any other 
technique in the science of materials. By measuring intensities of weak, space-group 
forbidden reflections, families of multipoles can be detected, including, magnetic charge 
(or magnetic monopole) dipole [5], anapole [El [7], quadrupole [8], octupole [91 flO] 
and hexadecapole pT| [T2] . Use of resonant x-ray Bragg diffraction to detect some of 
these strange multipoles hosted by complex materials is quite analogous to atomic parity 
violation experiments to measure the nuclear anapole [13] and it is the method by which 
to verify strange multipoles predicted in simulations [14] . 

A great deal of attention has focused on transition-metal materials especially 
manganites, in which it has been shown that orbital order is the key to understanding 
a raft of relevant physical properties [IHl [16]. The situation of rare-earth (4f) or 
actinide (5f) ions is different because the coupling between the spin and orbital angular 
momentum is relatively strong. High-order multipolar interaction may eventually 
manifest themselves in other subtle effects, such as lifting the degeneracy between single 
or multi-k structure, or independently of any magnetic order, the interactions between 
the 4f shells drive orderings in which their electronic density no longer respect the 
initial symmetries. In the latter cases, the primary order parameter is not of magnetic 
nature, but of orbital one. As the additional 4f asphericity is usually well described by 
the emergence of a quadrupolar arrangement, this modification of the system is called 
quadrupolar order. The case of multipolar ordering, although its possibility has long 
been recognized and theoretically investigated, remains a rather esoteric issue mainly 
because of the limited number of systems that unambiguously display the behaviour 
[1] . The situation has improved recently because novel experimental findings have been 
reported for several more rare-earth and actinides compounds [T7 | IT8 | |3]. 

The purpose of this article is to interpret resonant x-ray Bragg diffraction by UPds 
[T9| [20] . in particular the dependence of the intensities of the different super- lattice 
peaks on sample temperature and sample orientation (azimuthal angle) in terms of the 
possible lowerings of space-group symmetry and their associated quadrupolar structures. 
Russell- Saunders coupling scheme yields S = l, L = 5,J = L — 3 = 4. The large orbital 
moment gives rise to a strong coupling to the lattice and hence it is not surprising 
that this is a system in which quadrupolar effects are dominant. Because the U -U 
distance, which is around 4.11 A, is larger than that of the Hill limit (3.6 A), the system 
(UPds) behaves like a localised 4f-system. The low temperature specific heat of UPda 
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reveals a contribution from low-lying crystal-field (CF) levels, which gives rise to sharp 
phase transitions. The order would then results from a delicate balance between rather 
strong magnetic interactions and CF interactions. The possibility that quadrupolar 
order occurs in UPds was first raised in [21] and reiterated since then. 

McEwen and co-workers [22] re- analysed available experimental data for UPda, and 
produced both a new energy level scheme, with wave-functions, for the high-temperature 
phase, and a model for the progression of the phase transitions. 

2. Calculation of the RXS structure factor 

Resonant x-ray Bragg diffraction experiments [191 EQ] that we discuss were carried with 
the x-ray energy tuned to the Mjy edge of U at 3.728 keV. At this edge, the resonant 
scattering cross-section is dominated by electric dipole (virtual) transitions connecting 
3d3/2 and 5f5/2 electron states. The observation of satellite peaks at the M/y edge of U 
establishes both that long-range order of the 5f electrons occurs in UPds and that these 
experiments probe ordering of the quadrupolar moments in the 5f shell. 

For Bragg diffraction enhanced by an electric dipole event (El) the unit-cell 
structure factor, F, can be written as a product of a quantity X, which describes the 
condition of the x-rays, and a quantity ^' that is related to the valence electrons [23] 

F = 5^X^.v^^ = 5^(-1)«X5qV^$. (1) 

K K,Q 

Here, i^' = 0, 1, 2..., is the rank of a spherical component and the projection Q takes 
values within the range —K < Q < K. Xq is described in terms of coordinates fixed by 
the chosen geometry of the experiment, while \I/q is written with respect to axes in the 
crystal. \1/q , a structure factor of electrons in a unit cell, is a sum of multipoles {Tq)^, 
at sites d in the unit cell multiplied by the usual spatial phase factors, 

vl/$ = ^exp(zk-d)(T^)d. (2) 

d 

Point-group symmetry places restrictions on the allowed components of (Tq). 
Measured intensities reveal information on anisotropy in time-even, charge-like (K = 
even integer) and time-odd, magnetic (K = odd integer) distributions. For scattering 
enhanced by an El event one can observe multipoles with rank K = (monopole or 
charge), K = 1 (dipole) and K = 2 (quadrupoles). 

3. Properties of UPds 

At room temperature, the material crystallizes in the double hexagonal close-packed 
structure (dhcp) PGs/mmc with U ions at sites with locally quasi-cubic (2a Wyckoff 
positions) and hexagonal (2d Wyckoff positions) symmetry. The palladium ions are at 
6g and 6h positions. Benefiting from the selective character of RXS, Pd ions play no part 
and the observed RXS intensity is coming from quasi-cubic uranium ions, given the fact 
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that ordered quadrupole moments are predominantly at the quasi-cubic sites [211 IIH] • 
Resonant x-ray diffraction measurements [HI |20] were done at reffections (103) and 
(104) in ortho-hexagonal axes with lattice constants a = y/Sahex, b = ahex, and c = Chex- 
In order to align crystal coordinates with the x-ray coordinates defined in [23], it is 
necessary to perform a rotation of coordinates defined by the Euler angles (vr, /Sio;, tt) 
with /3io3 = 108.0 degrees and /3io4 = 103.7 degrees. Measurements in ^19j were only 
done at a fixed azimuthal angle tp= 90 degrees. 

Bulk measurements (heat capacity, electrical resistivity, ultra-sounds measure- 
ments, magnetic susceptibility, among others) together with x-ray and neutron scatter- 
ing experiments [21] have revealed up to four quadrupolar phase-transitions: Tq = 7.8 K, 
T+i = 6.9 K, T_i = 6.7 K and T2 = 4.4 K, being the transitions at Tq and T+i second 
order and the transitions at T_i and T2 first order. In the latter phase, the antiferro- 
quadrupolar (AFQ) transition is accompanied by the appearance of a small antiferro- 
magnetic dipolar moment. 



4. Minimal symmetry model structure-factor 

We start by deriving the consequences of a minimal symmetry model structure factor 
similar to the one used in refs. [25] and [26], and investigate the possibility of a 2-fold 
axis of rotation symmetry about the b crystallographic axis. The general structure 
factor in atomic axes rj, () with zeta parallel to the crystal b-axis, can be written as 

^Si^,V,0 = BS + D^ (3) 

with B^Q = Bq and B)^q = —Dq. Atomic axes {C,, are related to crystal 
b,c) by a rotation of tt/2 about the a-axis. If we calculate the resonant structure 
factor for rank 2 tensors (quadrupoles) in reflections of the kind (lOZ) at azimuthal angle 
of 90 degrees, we obtain, 

F.'. = ^JIb^o 

F^.„ = ReDlcos{(3 + 6)+ ImBlsin{(3 + 9) 

F^.^ = ReBl{2cos^l3 - 1) - -^Bl cos 26 + ImDl sin 2(5 

v6 

F,,^ = ReDl cos(/3 -6)+ ImBf sin(/3 - 9). (4) 

Assuming a point symmetry including a 2-fold axis about b-axis, \Q\ is an even integer. 
Thus Bf and Df would be forbidden leading to zero intensity in the rotated channels 
ttV and ct'tt at ip= 90 degrees. This is compatible with the RXS measurements done at 
7.1 K [201 [12] • In the measurements done at T= 5.2 K, the appearance of intensity in 
the rotated channel at the azimuthal angle ip = 90 degrees implies that the 2-fold axis 
point symmetry about b is no longer present. 



Antiferro-quadrupolar structures in UPd^ inferred from x-ray resonant Bragg diffractionh 



5. RXS measurements at T = 7.1 K 

There is evidence that the material undergoes a structural transition to an orthorhombic 
space group [271 [28]. Among the orthorhombic subgroups of PGs/mmc compatible with 
a 2-fold axis of point symmetry along b for uranium ions occupying quasi-cubic sites, 
we find that P222i reproduces the azimuthal measurements published in ref. [20]. In 
P222i the (10/) reflections would became space group allowed, in agreement with the 
weak non- resonant intensity in the unrotated channel reported in [19]. The quasi-cubic 
uranium ions would be at positions 2c and 2d of P222i. Uranium ions with z = 1/4 
and z = 3/4 are related by a two-fold axis about c. The quadrupolar structure factors 
for the 2c and 2d sites for the reflection (10/) and arbitrary Miller index / are, 

vl/|(10/)(2e) = exp(m/2/)(l + (-l)'+«))r2)(2,) 

^imm = - exp(m/2/)(l + (-l)'+«)(r<^)(2,). (5) 

In these structure factors, we use a different set of axes for the multipoles than those used 
in previous section, for the axes of (Tq) are now aligned with crystal axes a, b, c. At the 
reflection (10/) what it is measured is the difference between the quadrupolar moments 
in the (0,0,0) and (|, |,0) orthorhombic positions, namely, (Tq)(o,o,o) ~ i,o)- 
they occupy inequivalent Wyckoff positions (2c and 2d) quadrupoles at positions (0, 0, 0) 
and (|,|,0) are not related by the space-group. A quadrupolar ordering in which 
(Tq)(o,o,o) = 1 0) would not be observed at the (10/) reflection. 

The factor (1 -|- (—1)'+'^) in (5) implies that multipoles (Tq) with \Q\ = 1 are 
observed in the (103) reflection, while (104) would give access to multipoles with even 
\Q\ ( \Q\ = 0,2). In addition, point symmetry (with a 2-fold axis about b) forces (Tq) 
to be purely real, and this means that in the (103) reflection Iie{Tf) = (Qzx) is observed 
while at the (104) reflection Re(T2^) = (Qx^^yi) and < Tq >= {Q-iz^^r'^) are observed. 

For reflections of the kind (10/) with arbitrary / the total structure factor will be 
proportional to the difference of multipoles between the two crystal Wyckoff positions, 
(Tq-)(2c) — {Tq)(^2cI), i-e., to the difference of the mean values of the multipoles between 
the sites (0,0,0) and (|, |,0) in orthorhombic coordinates. 

The absence of measured resonant intensity at (104) at this temperature [TH] implies 
that (Qsz'^-r^) and (Qxz) have no antiferro-ordering. The intensity of (103) reflection 
is proportional to the difference of (Qzx) between the 2c and 2d crystal positions, 
{Qzx) {2c) - {Qzx) {2d) and is given by 

Fa'a{w?,) = - 2 {Qzx) sin (3q cos ao (6) 
^7r'<7(io3) = - (Q^x)[sin(ao + 7o-6')(l + cos/3o)(2cos/3o-l)) + 

sin(ao-7o + ^)(l -cos/3o)(2cos/3o + l))] (7) 

The angles ao, /3o; 7o result from the composition the azimuthal rotation and the rotation 
that aligns the Bragg wavevector with the x-ray coordinates, and are given by 

ao = arccot (— cot ip sin (3) 
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Figure 1. Azimuthal variation of tlic intensity in the (103) reflection measured at T= 
7.1 K. Experimental points are taken from ref. 20 . The continuous line shows the fit 
of the data to eqs. (6) and (7). 



/So = arccos (cosV' COS/3) 

7o = arccot (— cot (3 sin ip) (8) 

with the angle (3 for the (103) and (104) reflections being: /3io3 = 108.0 degrees, 
/3i04 = 103.7 degrees. 

The azimuthal dependence of the (103) reflection together with experimental 
data [20] is depicted in Fig. 1. We have allowed contribution from a non- resonant 
(Thomson) term in the a'a channel. The quadrupolar structure depicted in fig. 4 
of ref. [20] corresponds to the antiferro-quadrupolar ordering of {Qxz) observed, i.e. 
{Qxz){2c) = —{Qxz) (2d)- A possible contribution of {Qx^-y^) that would modify the 
azimuthal dependence of the (103) reflection was proposed in [29]. However, in space 
group P222i contribution of (Qx^-y^) to the (103) reflection is not allowed, as it would 
have in-phase stacking along the c-axis, i.e. when making the translation (0, 0, |) along 
z, the value of {Qx'^-y^) would be the same (this is different to the structure depicted 
in [12] which has antiphase stacking of {Qx2^y2) along c-axis). Still, in the structure of 
space group P222i there would the possibihty of ordering of (Qx^-y^) with {Qx^-y'2){ooo) 
= {Qx^-y^) (L^ifi), which would not be observable in (10/) reflections. 
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6. RXS measurements at T = 5.2 K 

To accommodate the observed change to the shape of the (103) azimuthal curve, and the 
appearance of intensity in vrV at 90 degrees in the azimuthal curves of both (103) and 
(104) reflections (see Figs. 2 and 3), it is necessary to reduce the symmetry experienced 
by the quasi-cubic uranium ions (two-fold axis about b). Possible subgroups of P222i 
that do no break the translational symmetry are P2i, P2, PI. The possibility of a 
lowering of symmetry to a monoclinic space group has already been suggested in ref . |28] . 
The minimum reduction of symmetry to justify these changes would be to consider the 
monoclinic space group P2i. In this space group, ions at different positions in z within 
the unit cell continue to be related by a two-fold axis about c. The quasi-cubic uranium 
ions in this group would occupy two different sets of 2a Wyckoff positions. The algebraic 
form of unit-cell structure factors is unchanged, but now there is no point symmetry 
and (Tq) is no longer purely real. In consequence, there are additional contributions 
containing (Qyz) at (103) and {Qxy) at (104). The azimuthal dependence of the (103) 
reflection is then: 

^a'<7{io3) = - 2sin/3o(((5^x) cosao + {Qyz) sinao) (9) 
^7rV(io3) = - {Qzx)[sm{ao + 'jo -0){1 + cos/3o)(2cos/3o - 1)) + 
sin(ao - 70 + ^)(1 - cos/3o)(2cos/3o + 1))] + 
{Q yz) [ cos{ao + 'Jo - 6')(1 + cos/3o)(2cos/3o - 1)) - 
cos(ao - 70 + ^)(1 - cos/3o)(2cos/3o + 1))], (10) 
and for the (104) reflection, 

Fa'aiW'i) = {Q zz) {3 COS^ f3o - 1) 

+ 2 sin^ f3o{^{{Q XX - Qyy)) COS 2ao+ (Qxy) sin 2ao) (11) 
3 

FTr'aiio4) = - -jkQzz) slu 2/3o sin(7o - 6*) 
sin /3o(l + cos /3o) 

{^((Qxx - Qyy))W + cos/3o) sin(2ao + 7o - ^) + 
(1 - cos /3o) sin(2a;o - 7o + Q)\ - 
(Qxj/)[(1 + cos/3o) cos(2ao + 7o - 6*) + 

(l-cos/3o)cos(2ao-7o + ^)]}. (12) 

Fig. 2 shows the fitting of the available data for (103) [20] with (Qzx) and (Qyz) in (9) 
and (10) allowed to be different from zero, together with the other quadrupoles. The 
value of {Qzx) / {Qyz) from the fitting is: (Qzx) / {Qyz) = —4.72 ± 0.12. Measurements 
done with vr incident polarization [19] show that at T=5.2 K, at tp=90 degrees, in (103) 
tt'tt is about one order of magnitude more intense than a'lr (with our derived parameters 
we obtain 1^,^/1^,^ = 19.4 ai ip = 90 deg. and I-k'-k / la'-K = 8.3 at ip = —90 deg.). 
This differs from the fit done in [20], which attributed the variation in the azimuthal 
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Figure 2. Azimuthal variation of tlie intensity in the (103) reflection measured at T= 
5.2 K. Experimental points are taken from ref. [20]. Continuous lines show the fit of 
the data to eqs. (9)-(10). 

dependence in (103) at 5.2 K solely to the contributions of {Qxy) and {Qx'^-y^)- In our 
proposed structure these two quadrupoles have in-phase stacking along c-axis and do 
not contribute to (103) reflection. 

Available data [201 EO] for the (104) reflection is shown in Figure 3 together 
with the azimuthal dependence calculated from (11) and (12) with only {Qxy) allowed 
to be different from zero and a non-resonant contribution (Thomson) in a'a that is 
independent of the azimuthal angle. Inclusion of (Qx^-y^) and (Qzz) does not lead to 
a significant improvement in the fit to (104) data. Moreover, in the case of the (104) 
reflection measurements done with vr incident polarization at T = 5.2 K [19] show a 
very weak intensity in the tt'tt channel and both (Qx^-y^) and (Qzz) give a non-zero 
intensity at 1^=90 degrees in tc'tt. However, an accidental cancellation could occur if 
XX Qyy) '^{Q zz} ■ 

Figure 4 shows the antiferro-quadrupolar orderings of (Qyz) and {Qxy) derived from 
the proposed space-group structure, which are compatible with data collected for the 
(103) and (104) reflections at T = 5.2 K and coexist with the ordering of (Qzx) observed 
at 7.1 K [20j. The stacking along the c-axis of the (Qyz) component of the quadrupolar 
moment shown in Fig. 4-a changes sign when making a translation of 1/2 of the z 
coordinate, which is different to the (Qyz) stacking used in refs. [20] and [30] . 
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Figure 3. Azimuthal variation of tlie intensity in the (104) reflection measured at T= 
5.2 K and 2.1 K. Experimental points are taken from ref. [30]. Continuous lines show 
the fit of the data to eqs. (11) and (12). 



4-a) Oi-dering of Qxy observed 4-b) Ordering of Qyz observed 

at T=5.2 K in (104) reflection T=5.2 K in (103) reflection 




Figure 4. Antiferroquadrupolar orderings of {Qxy) and {Qyz) observed in the (104) 
and (103) reflections at T = 5.2 K. 



Antiferro-quadrupolar structures in UPd^ inferred from x-ray resonant Bragg diffractionlO 

7. Conclusions 

Resonant x-ray Bragg diffraction data are evidence of an antiferro-quadrupolar ordering 
in UPda- Our analysis of available data identifies the space group for the phase between 
To = 7.8 K and T+i = 6.9 K as P222i. Additionally, we can infer the quadrupole 
components (Qa/s) («, (3 — x,y, z) and their stacking pattern contributing to different 
space-group forbidden reflections. Intensity in the (103) reflection is here assigned to 
the (Qzx) quadrupole. For the phase below T_i = 6.7 K, measurements are compatible 
with the space group P2i, which implies loss of the 2- fold axis of rotation symmetry 
in the point symmetry sites used by the U ions. However, we should be cautious with 
this assignment of space group, given that the phase transition at T_i is first-order and 
in that case the low temperature space-group is not necessarily a subgroup of P222i. 
The intensity at the (103) reflection would come mainly from (Qzx) together with a 
small (Qyz) contribution. Intensity at (104) is attributed to {Qxy) with a possible small 
contributions from {Qx2_y2) and (Qzz)- Within the derived structures for both phases 
of the compound, no additional contributions from other quadrupolar parameters are 
present in (103) and (104) reflections. 
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